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The reduced dynamics of an atomic qubit coupled both to its own quantized center of mass 
motion through the spatial mode functions of the electromagnetic field, as well as the vacuum 
modes, is calculated in the influence functional formalism. The formalism chosen can describe the 
entangled non-Markovian evolution of the system with a full account of the coherent back-action of 
the environment on the qubit. We find a slight increase in the decoherence due to the quantized 
center of mass motion and give a condition on the mass and qubit resonant frequency for which the 
effect is important. In optically resonant alkali atom systems we find the effect is negligibly small. 
The framework presented here can nevertheless be used for general considerations of the coherent 
evolution of qubits in moving atoms in an electromagnetic field. 
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Atomic motion is an unavoidable element in the consideration of any AMO system and an integral part of 
experimental designs in atom trapping devices. At issue here is the interaction between the internal degrees of 
freedom of an atom, assumed to contain an effective two- level system (qubit), and the electromagnetic field (EMF), 
modified by the atom's quantal motional degree of freedom. This problem has two aspects: 1) How does the two level 
activity affect the atomic motion? and 2) How does atomic motion affect the two level activity? The first aspect is 
the basis for laser cooling and atom trapping, which have been studied in great detail and successfully implemented 
by well-known experiments (for reviews see pi 0, 0, 3 ) • This paper is aimed at the second aspect, specifically, how 
quantized motion affects the qubit-EMF system dynamics, which is of interest in the design of quantum computers 
based on atomic qubits (in the form of a neutral atom 0, El 0> H> 13 or i° n US) in a QED cavity or optical potential. 
Effects on internal dynamics due to quantized center of ma ss ( COM) motion have previously been studied in the 
situations of an atom in free space > m a cavity [l2l Il3l Il4| , and when the atom's qubit and COM degrees of 
freedom are entangled |l5j . However, all have focused on spontaneous emission rather than decoherence. The present 
O |' work probes the non-Markovian regime of atom-EMF interaction, under the modest aim of explicitly computing how 
entanglement with quantized motion through recoil affects the decoherence and relaxation rates of an atomic qubit 
in free space. In order to achieve that end, we first discuss two issues of importance in computations of coherent 
Ci \ reduced dynamics, using path integral methods. 
2 . 

i The importance of including back-action 

It is well-known that the interaction between a two-level system (2LS, or qubit) and the EMF is the primary 
source of its relaxation and decoherence, while effects associated with the atom's motional degrees of freedom are 
usually relegated to the background. Assuming that the atom moves adiabatically limits one's consideration to those 
circumstances wherein the external degrees of freedom act merely as a passive parameter in the environment (here 
comprised of the EMF and atomic motion) of our system (the qubit), with no dynamical interplay. In technical 
terms, this amounts to a 'test-field' approximation - that the qubit lives in a fixed environment defined by a set of 
parameters, amongst them the adiabatic motion 27] . The test field approximation leaves out effects of changes in 
the environment on the system. To include the effects of the environmental variables dynamically it is essential to 
perform a self-consistent back-action calculation. This was done for the effect of a cavity EMF on the 2LA in Ref. [l6| . 

Full coherence requires self-consistent treatment 

In tackling problems where many factors enter, it is useful to isolate one factor after another so that the re- 
maining factors of interest to us can be simplified enough to yield some solution. For quantum coherence and 
entanglement such simplifications can lead to erroneous results, since phase information is lost if one artificially 
isolates the linking components of the complete quantum system. This brings up the necessity of self-consistency 
in any treatment of quantum coherence and entanglement issues. In the present case of a qubit in an EMF this 
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requires that the fully entangled system of atomic 2LS (internal), the EMF, and the center of mass (external) 
degrees of freedom be treated coherently as a whole and each factor involved be allowed to evolve under the 
influence of the others in a self-consistent manner. This self-consistency requirement leads to non-Markovian dy- 
namics since memory effects arise naturally and are necessary to preserve maximal coherence during the evolution [28j . 

Grassmannian and Coherent State Representation of Influence Functional 

A theoretical scheme we found satisfactory in meeting these requirements is the influence functional (IF) formalism 
of Feynman and Vernon or the- related closed-time-path (CTP) effective action of Schwinger and Keldysh [l8| . 
The influence of the environmental variables on the system of interest is incorporated in the IF (or effective action) 
in such a way that the equations of motion obtained for the system will already have included the back-action of 
the environmental variables on the system in a dynamically self-consistent manner. This scheme has been applied 
to a two-level atom (2LA) interacting with an electromagnetic field (EMF) in reference 0]. There, a first-principles 
derivation of the general master equations is given and applied to the study of the decoherence of a 2LA in an EMF, 
for the cases of a free quantum field and a cavity field in the vacuum at zero temperature. The authors found that for 
the standard resonant type of coupling characteristic of such systems the decoherence time is close to the relaxation 
time. 

Here we use the influence functional method for the treatment of the back-action of the quantum field and the 
quantal motion of the atom on the qubit. In Section II we compute the transition amplitude between an initial and 
final state using a coherent state label for the (bosonic) states of the EMF and a Grassmannian for the (fcrmionic) 
2LS. The coherent state basis allows us to identify the Hilbert space of states with a space of coherent states. The sum 
over all quantum evolutions is then a sum over all paths in this space. Once the transition amplitude is computed in 
some sufficiently simplified form, forward and backward versions can be combined and reduced to form the reduced 
density matrix evolutionary operator. In Section III we calculate the evolutionary operator for the reduced density 
matrix when the EMF and motional degrees of freedom are integrated over. We derive an equation describing the 
evolution of the on and off-diagonal elements, the latter is the coherence function we seek. We end in Section IV with 
a discussion of our results and comments on possible further developments on this subject. 



II. THE TRANSITION AMPLITUDE 



Our system is a 2-level atom interacting with its own center of mass (COM) motion and the EMF. We begin with a 
modified multi-mode Jaynes-Cummings type Hamiltonian (e.g. dipole and rotating wave approximation, see Appendix 
A of Ref|T3), 

p2 

H = ^j + h ^S z + hJ2i^blb k + g k (X)S+b k + g k (X)S_4]. (1) 

k 

The first term in the Hamiltonian is the COM kinetic energy. The next two terms are the qubit and EMF energies, 
respectively. The last two terms are the interaction between the qubit, EMF, and the atom's COM degree of freedom. 
Note that P and X are both operators. Coupling of the qubit to its COM motion is through the spatial mode 
functions of the EMF. We shall restrict our consideration to an initial vacuum EMF at zero temperature. The result 
of this calculation will thus be the modification of the vacuum decoherence and relaxation rates of a qubit when the 
effects of quantized atomic motion are included. 

The first step towards obtaining the reduced system dynamics while retaining the full system's coherence is to 
compute the transition amplitudes between the initial and final states which are the matrix elements of the evolution 
operator of the full system. We do this with coherent state path integrals. For the EMF we use a bosonic coherent 
state representation and for the 2-level system (qubit) degree of freedom we use the Grassmannian coherent states 
[20I l2l) . Coherent states are by definition generated by the exponentiated operation of the creation operator and a 
suitable label on a chosen fiducial state: 

|z k ) = exp(z k 6 k )|0 k ) (2) 
\rj) = eMvS+)\0) (3) 

In the case of bosonic coherent states defined in Eq. © the label, z k , is a complex number, and in the case of the 
Grassmann coherent states defined in Eq. the label, 77, is an anti-commuting number. The chosen fiducial states 
are the EMF vacuum and the lower 2-level state, respectively. 

In order for any set of states to be useful for the decomposition of the transition matrix they must have a resolution 
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of unity. The EMF and Grassmannian coherent states have the following decompositions of unity 

i = / d M (* k )kk><2k| = / Mv)\v){n\ (4) 



with the measures 

dfj,(zb) = exp(-z k z k ) 
dn(rj) — exp(—fjr)) 

Grassmann coherent states also share other well known properties of coherent states such as being non-orthogonal 
and eigenstates of the annihilator: 

(zk|z' k ) = exp(z k z' k ) {fj\rf) = exp(fjr)') 
bk\z k ) = z k |z k ) S-\rj) = rj\r}) 

The center of mass or external degree of freedom can be represented in either the position or momentum basis. In 
the coherent state basis the Hamiltonian Eq. (JIJ can be written in its Q-representation [22|, |23|, yj, |2f| as [cf Eq. (2.8) 
of [3] 

^({^k},{2k}, j?,?7,X) = — \- hu fjrj + Ti >^[a;k^kgk + ?7gk(X)0k + ^kgk(X)7?]. (5) 

k 

The transition matrix elements between the initial and final coherent states are then 

K(t,0) = ({^kj^X^tlexp^i^Kz^X^O). (6) 

Using the completeness property of the (EMF and Grassmann) coherent state basis to facilitate time-discretization 
of the transition matrix 17] puts the transition matrix elements in a coherent state path integral representation. After 
inserting the Q-representation, the transition elements transform into a sum over paths in the coherent state labels. 
Having done the above the transition matrix becomes a triple functional integral: 



K(t,0) = J DltJ DfjDr,Y[Dz k D Z]i eMfifV(t) + ^Z/kztf) J X da] e 1 



iui t/2 



x exp 



/ (7777 + iw 7?77 + ^z k z k + i^w k z k z k i^?7g k (X)z k -H^z k g k (X)?7)ds 

1, 1, 1, 1, 



(7) 



In this form the transition matrix elements can be evaluated exactly by a combination of stationary phase and 
correlation function methods which exploit the truncating properties of Grassmann variables. The order of evaluation 
will be the EMF, COM, and then Grassmann functional integrals. The details follow. 

A. EMF Path Integral 

First, the EMF coherent state part of the triple path integral can be evaluated by the stationary phase method 
(I?! . The variational equations of motion for the electromagnetic field variables in Eq. Q are 

z k = -iw k z k - ig k (X)?7 (8) 

which have integral solutions [cf Eq. (2.14) of 16]] 

z k (s) = z lW e'^ s - i f dr g k (X(r))e-^ s - r ^(r). (9) 



(i 



The transition amplitude from an initial EMF vacuum ({zi k } = 0) to an arbitrary final state becomes 
K(t,0) = J DfjDrjexp fj f r)(t) - J (fjr) + iw fjr) + ^X )ds 

x JJexp i 1 ns- Ak [./\.\si .]< . /k 

k 

The path integral for the EMF degrees of freedom is now complete 



e iWot/2 

i / dsg k (X( S ))e- i ^ t - s ^ /k ry(s) - f da f drg k (X( S ))g k (X(r)) e -^( s - r )^( S )r / (r) 

Jo 



.(10) 
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B. COM Path Integral 



Second, the position path integral can be evaluated as a set of 0, 1 and 2 point functions. Note in the transition 
amplitude of Eq. i|l(Jfl that since the EMF is taken to be in an initial vacuum, any source term for rj(s) will be 
proportional to rji. The variational equation of motion derived from Eq. I|1U|) for ?y(s) yields 



)(s) = -iUoV(s) ~ i / dr Vg k (X(s))g k (X(r))e" 
Jo 



rj(r) 



(11) 



with the boundary condition 77(0) = r\i. Therefore 77(5) = u(s)r)i. We use this to expand the exponent in the transition 
amplitude of Eq. (|10Jl . Due to the nilpotency of the Grassmann variables (i.e. rjf = 0) it will truncate after the first 
term in the expansion. 

After expanding and truncating the integrand, the position path integral is 



/ 



DXexp 



(12) 



1 - i / t d S ^g k (X( S ))e-^( t - s )z /k r,( S )- fds / S dr^g k (X( S ))g k (X(r))e- 
Jo ^ Jo Jo v 



There are thus three correlation functions which need to be computed. First the spatial mode functions must be 
chosen in order to specify the targeted correlation functions. For an electromagnetic held in free space (no cavity or 
boundaries) 



gk(X) 



A 



/LUk 



: exp(ik • X) 



(13) 



the correlations functions are computed in Appendix A. Substituting these expressions back into the Eq. (|1(J|I gives 
for the transition amplitude 



K(t,0) = I DfjDTjexp 



lj\£ 

— (X f -X,) 2 +fj f rf(t) - / (fjfi + iu fjri)ds 



A 



exp 



i-k - (X, 



(—) 

\2irit J 
i s(t-s) 2 



3/2 



e iw t/2 e -iw k (s-r) 



2M t 



rt rs ^2 

/ ds dr) — exp 

J0 JO V w k 



i- — -k- (X, -Xi) 

t w ; 2M 



iuJk(t - s) z fk rj(s) (14) 
fj{s)r)(r) 



i (t-( s - r ))( S -r), 2 



The path integral for the external degrees of freedom is now complete. 



C. Qubit Path Integral 



Finally, the Grassmann variable path integral can be evaluated along its stationary path. The variational equation 
of motion for the Grassmann held variable in Eq. I|14(l is 



f s A 2 
Vt(s) = -U0 o rj t (s) - / dr — /j t 
Jo V w k 



(s - r) r) t (r) 



(15) 



with the definition: 



Ht(s) = exp 



iwkS — i-k • (Xf — X,) 



i S(t-S) k2 



(16) 



Note that the hnal time t enters as a parameter in the variational equation of motion just as the mass or position do. 
The reason for this is that the above variational equation of motion is for the evolution of the atom from an initial 
time to a final time, so the time is an explicit parameter. 
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Rewriting the above variational equation in Laplace space allows the non-local integral part to be transformed with 
the convolution theorem. The solution is in terms of an inverse Laplace transform, 



Vt(s) = ViUt(s) 



Th_ n +io ° e sz dz 
27ri J 7 _ ioo z + \uj + p(z) 



with the definition: 



fit(z) = — / e sz exp 

Jo 



i s(t-s) k2 



dz. 



(17) 



(18) 



The solution thus becomes a contour integral. The pole of the denominator in Eq. (|17|l can be found to 0(X 2 ) 

Z o = -l0J o - /i(-iw ) +0(A 4 ). (19) 
Finding the pole to order 0{\ 2 ) gives a solution to the same order: 



u t {s) = e" 



exp i — A 2 ty^— / dsexp 
I 4^ J 



1 \0Ju - U> + 



k- (X 



t 



2M 



ik 2 

2 Aft 1 ' 



(20) 



Evaluating the transition amplitude along its stationary path with the second order pole approximation yields an 
expression for the transition matrix that is second order in its action: 



/ M \ 3/2 

A '«-°>=M exp 



i M 

iL0 o t/2 + — (X / -X,) 2 + 0(A 4 ) 



(21) 



exp 



VfVt 



(t) -if ds —— exp [ - iw k 

JO ^ V W k 



-iw k (t-s) + i-k-(X / -X i 



is(t - s) 
2Mt 



k 2 ]z/ k r?t(s) 



All three functional integrals are now evaluated. In the next section we proceed to derive the evolutionary operator 
for the density matrix by combining the transition amplitudes into a closed loop. 

III. EVOLUTIONARY OPERATOR 

At this point the expression of Eq. (|21|) for the transition amplitude can be combined with its counterpart prop- 
agating backwards in time and traced over all final EMF states. The result gives the evolutionary operator for the 
reduced density matrix (we may call it the reduced propagator, for short), 



J R = J dX f ]Jdn(z f k)K(t,0)K*(t,0), 



(22) 



and is formed by integrating out the environmental variables which in our case are the EMF and the atom's motional 
degrees of freedom. 

The evolution of the qubit density matrix elements with back-action from the EMF and the atomic motion can be 
calculated from the reduced propagator 

Pr® = J d^d^dfiiX,) J R p A (p) ® p x (0). (23) 

The functions pa(0) and /Ox(0) are initial states for the 2-level atomic and external degrees of freedom, respectively: 

p A (t) = poo(0) + w>io(0) + »?-poi(0) + i^pii(O) (24) 



P xW- < f > (x i )o>*(x l ) 



(25) 



The function $(X) is the initial (external) center of mass wavefunction of the atom. From Eq. ll'ol) the on and 
off-diagonal components of the reduced density matrix elements evolved to time t are given by 



Pn{t)= Pn (0) 



M 



-Y 

2?rit / 



dX, / ax; / dX l $(X l )$*(X^)exp|^(X / -X l ) 2 -^(X / 



u t {t)u t {t) (26) 



6 



0.20 P 



° - \ 

CD - \ 

u 0.10; 
c . \ 

K ° 0.05 ; \ ; 

0.00 I T~r~ - 

3 4 5 6 7 8 
log 10 m* 

FIG. 1: A plot of the percentage increase in the decoherence of the off-diagonal matrix elements of the reduced density matrix 
versus the non-dimensionalized mass (to* = tt 2 -)- The decoherence rate increases as the mass of the atom containing the 
qubit is decreased. As the mass is increased the decoherence rate asymptotes to the value of a stationary atom obtained by 
Anastopoulos and Hu [Tfij |. This is consistent with a smaller mass qubit being more affected by its recoil than a heavy mass 
qubit. Typical experimental parameters fall to the right end of the shown plot. 



PloO) = Pio(O) 



JL) 



J dX, J dX^ J dX. i $(X. i )$*(X' J )exp|^(X / 



iM 



(X/-XJ)' u t (t). (27) 



The EMF, as previously stated, is in a vacuum state, but the choice of an initial center of mass wavefunction has 
not yet been made. To closely model an atom with fixed position and momentum, we use a minimum uncertainty 
Gaussian wavefunction centered at (X Q = 0, P Q = 0). 



*(X) = ^- 3 / 4 a- 3 / 2 exp 



2d 2 



(28) 



Such an initial wavefunction simplifies the expressions for the diagonal and off-diagonal matrix elements of the qubit. 

The result for the off-diagonal components which measures the coherence of the qubit under such conditions is 
shown here: 



PioW = Pio(0) 



4 / M 2 a 2 \ 3/2 f 00 , 2 , . 



M 2 a 



2^-2 



t 2 - 2iMa 2 f 



(29) 



The function u(x,t) is given by Eq. I|20l) with x = |X/ — X»|. 

The evolution of the coherence function is found to follow an exponential decay with a decay rate slightly faster than 
in the infinite mass case. The percentage change in the decoherence rate of the off-diagonal versus the the stationary 
qubit case is plotted in Fig. (1). The decay rate increases with decreasing mass and matches the stationary qubit 
result given by m the limit of infinite mass. We expect that a qubit in a smaller mass object is more affected by 
recoil than a qubit in heavy mass. The variation in the decoherence rate with changes in the external wavefunction 
size is relatively flat and cannot reliably be resolved with the available computing power and machine accuracy. We 
find that so long as the resonant frequency is small enough or the mass large enough that the atomic recoil velocity 
is non-relativistic, which is where this theory is valid, then the motional decoherence will contribute negligibly to the 
decay of the qubit. 



IV. DISCUSSION 



Often, one may separate the dynamics of an atom's motion from those of its internal degrees of freedom by arguing 
that the time scales associated with the motion of the atom are much longer than those of the two level activity. This is 
the rationale behind the adiabatic approximation adopted for most considerations of the atomic dynamics. However, 
coherence requirements in quantum computing implementations may prompt one to question this assumption. One 
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aim of our investigation is to test for non-adiabatic effects in atomic quantum computing schemes. Another is to 
describe the effect of recoil from the emission and re-absorption of virtual particles in the atom-EMF interaction upon 
the center of mass motion. These two problems correspond to the two aspects described in the Introduction. Here 
we consider the second aspect mentioned above, aiming at the effect of quantum motional decoherence of the qubit, 
i.e., the back-action of atomic motion on a two level system in free space as mediated by the EMF. 

We find that the inclusion of the external degrees of freedom only slightly alters the decoherence and relaxation 
rates as compared to a stationary atom. Typical experimental parameters fall to the right end in the plot of Fig. (1). 
A Rubidium atom used as a qubit would have a non-dimensionalized mass of approximately Log(-^j-) = 8, which 
places it in a regime in which the effect of motion-induced decoherence is negligible. For optical qubit transition 
frequencies in general, motion-induced decoherence will not be a factor unless the mass of the qubit is four to five 
orders of magnitude smaller than the mass of a typical alkali atom. One can conclude tentatively that in general 
AMO implementations, motion-induced decoherence of a free qubit is negligibly small. Since the calculation done 
here is coherent and non-Markovian, one can view our result as confirming the validity of the adiabatic approximation 
in alkali atom qubits. 

Although the result of the calculation is the expected one, the technique described here is the first able to compute 
the decoherence of a qubit coupled to its own quantized COM without any form of Markovian approximation, while 
allowing the qubit-EMF coupling to be the non-linear form derived from the EMF spatial mode functions. Useful 
applications of this method will include any situations in which the COM motion of an atom back-acts onto its 
internal qubit dynamics and the full multi-mode structure of the EMF is relevant. Two such examples, as drawn 
from the references cited in the Introduction, are an atomic qubit in a cavity and an atom with entangled qubit 
and EMF deg rees of freedom. In the former, the presence of the cavity walls increases the cavity mode recoils on 
the atom |2fij |. The latter is at the center of certain two qubit gate implementations Q, with the question there 
being how well coherence is maintained when a qubit is entangled both internally and externally. Calculation in 
that case can provide an important feasibility test of quantum computing applications which utilize such entanglement. 
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APPENDIX A: COM FUNCTIONAL INTEGRAL 



The position path integral which needs to be evaluated is: 



DXexp 



x exp 



iM 



ds X 



(Al) 



if d S ^g k (X( S ))e-^(*- s )% k r ? ( S )- / ds rdr^g k (X( S ))g k (X(r))e- i 
Jo Jo Jo ,, 



ia>k(s— r) 



which can be expanded and truncated to: 



/ 



DXexp 



iM 



ds X 



(A2) 



I - i rds^g k (X(s)) e ^^ s )z /k rKs)- ['da f dr £g k (X( S ))g k (X(r)> 
Jo Jo Jo ,, 



-ia> k (s— r) 



There are thus three correlation functions which need to be computed. First the spatial mode functions must be 
chosen in order to specify the targeted correlation functions. For an electromagnetic field in free space (no cavity or 
boundaries) 



A 



the correlations functions are: 

rX/.t 
Xi,0 



gk(X) 



DXexp[-^ / drX(r)] 



exp(ik • X), 



(—) 

\2TTlt I 



3/2 



exp 



iM nc 



(A3) 



(A4) 



i-Xf.t 



/x,,o 



DXexp[ik.X( S )-^rd.X ( .)]=(^) 



3/2 



exp 



(X, - Xi ) 2 + j£ k- ( X/ - X,) - 4- s(t ~ s) ) 



(A5) 
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r X/,t 



iM 



DXexp[-ik- X(s) + ik-X(r) 
/Xi,o 2 Jo 

/ M \ 3/2 



\27rit / 



exp 



dr X(r 
iM 



(A6) 



2t 



i (<-(g-r))(a-r) 
2M i 



Substituting these expressions back into the Eq. (|10J) gives for the transition amplitude: 



K(t,0) = / DfjDijexp 



+ Vflit)- I (r}fi + vu) a r)r))ds 
A 



(—) 



3/2 



exp 



iM 



1 - i 



Jo 



exp 



i S(t-S),_ 2 



i-k- (Xf-Xi) 
t y 1 J 2M t 



2t 



—iuj k (s-r) 



t r-s y2 

ds dr) — exp 
o Jo k Wk 



i- — -k -(Xf-Xi) 

t v 1 ' 2M 



(X/-X,) 2 
i (t — (s — r))(s — r) 



t 



(A7) 
fj{s)r](r) 
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